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^—\ , Context. Active galactic nuclei show a wealth of interesting physical processes, some of which are poorly understood. In a broader 

r~^ . context, they play an important role in processes that are far beyond their immediate surroundings, owing to the high emitted power. 

f^ ' Aims. We want to address a number of open questions, including the location and physics of the outflow from AGN, the nature of 

the continuum emission, the geometry and physical state of the X-ray broad emission line region, the Fe-K line complex, the metal 
abundances of the nucleus and finally the interstellar medium of our own Galaxy as seen through the signatures it imprints on the 
X-ray and UV spectra of AGN. 

Methods. We study one of the best targets for these aims, the Seyfert 1 galaxy Mrk 509 with a multiwavelength campaign using 
k> ; five satellites (XMM-Newton, INTEGRAL, Chandra, HST and Swift) and two ground-based facilities (WHT and PAIRITEL). Our 

rN , observations cover more than five decades in frequency, from 2 /jm to 200 keV. The combination of high-resolution spectroscopy 

^ ■ and time variability allows us to disentangle and study the different components. Our campaign covers 100 days from September to 

December 2009, and is centred on a simultaneous set of deep XMM-Newton and INTEGRAL observations with regular time intervals, 
spanning seven weeks. 

Results. We obtain a continuous light curve in the X-ray and UV band, showing a strong, up to 60% flux increase in the soft X-ray 
band during the three weeks in the middle of our deepest monitoring campaign, and which is correlated with an enhancement of the 
UV flux. This allows us to study the time evolution of the continuum and the outflow. By stacking the observations, we have also 
obtained one of the best X-ray and UV spectra of a Seyfert galaxy ever obtained. In this paper we also study the effects of the spectral 
energy distribution (SED) that we obtained on the photo-ionisation equilibrium. Thanks to our broad-band coverage, uncertainties on 
the SED do not strongly affect the determination of this equilibrium. 

Conclusions. Here we present our very successful campaign and in a series of subsequent papers we will elaborate on different aspects 
of our study. 

Key words. Galaxies; active - quasars: absorption lines - X-rays: general 



1. Introduction 

Active galactic nuclei (AGN) are compact sources with very 
high luminosities, located at the centres of galaxies. Accretion 
onto the super-massive black holes (SMBH) at their centres is 
generally believed to be the driving process for the activity. 
Thanks to their brightness, they form one of the richest labo- 
ratories for studying astrophysical processes. In this paper we 
present one of the deepest multiwavelength campaigns of an 
AGN, the Seyfert 1 galaxy Mrk 509. It is among the best for 
these studies because it is unique in combining X-ray bright- 
ness, outflow features, and significant but moderate variability. 
Below we introduce the most important astrophysical processes 
that are addressed by our study. 

After our introduction of the relevant astrophysics, we briefly 
provide an overview in Sect. [3] of the target of our campaign, 
Mrk 509, followed by an overview of the observations in Sect.lH 
and we present the spectral energy distribution (SED) and light 
curve in the subsequent sections. Details about various aspects 
of our study are deferred to subsequent papers of this series. 

2. Astrophysics of AGN in the context of our 
campaign 

2.1. Outflows from AGN 



In recent years 
flows for the 
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2001), and the luminosity function of AGN (iWvithe & Loeb 
2003h has been widely recognised. However, for the lack 
of a better alternative, theoretical studies use the physical 
properties of the outflow (metallicity, mass, and kinetic energy 
flux) as free parameters because there are few observational 
constraints. To assess the importance of AGN outflows on the 
processes mentioned above, it is essential to establish the real 
mass flux, metallicity, and kinetic luminosity (E^) of AGN 
winds. Determining ZJk requires answers to some fundamental 
questions. What is the physical state of the outflowing gas and 
what is its total column density? What is the distance of the gas 
from the central source? 

The iJk of a shell-like, non-accelerating outflow is given by 

Ek - ^nf^RNnmy, (1) 

where Q. is the solid angle occupied by the outflow, R the dis- 
tance from the central source, A^h the total hydrogen column den- 
sity, nip the proton mass, /,„ = p/nn'Wp (1.43 for a plasma with 
proto-solar abundances) with p the mass density and tin the hy- 
drogen density, and v the outflow velocity. Spectral observations 
straightforwardly determine v, and Q. is expecte d to be ~ ;r since 
50% of all Seyfert Is show outflow signatures dCrenshaw et alJ 
11999 ). Our campaign focuses on determining the two most un- 
certain quantities: A^h and R. Determining these is tightly con- 
nected to the physical state and location of the outflow, as ex- 
plained below. 

In the best studied AGN at least two to three ionisation com- 
ponents are needed to model the rich X-ray absorption spectrum. 



In NGC 3783 there are three discrete components, which are 
most likely in pressure equilibrium, representing different phases 
of gas at the same dista nce from the nucleus ( Krongold et al] 
i2003L iNetzer et al.ll2003h . However, in NGC 5548 at 'least five 
ionisation components are needed if the X-ray absorber is mod- 
elled by a finite number of discrete absorption systems, and 
these absorbers cannot be in pressure equilibrium ( Kaastra et al.l 
2002). Instead, a continuous, power-law distribution of A^h ver- 
sus ionisation parameter f give s a better description of the data 
feteenbrugge et al.ll2003T l2005h . In yet another case, Mrk 279, 
the distribution is continuous but mor e complicated than a sim- 
ple power-law (C ostantini et al.ll2007l) . 

Photo-ionisation modelling of the outflow yields the ionisa- 
tion parameter ^ = L/nR^, where L is the 1-1000 Ryd ionis- 
ing luminosity and n the hydrogen density. Independent mea- 
surements of n then determine R. The use of density sensitive 
X-ray lines is difficult and has not (yet) delivered robust results 
(iKaastra et al.l2004.) . Alternatively, when L changes, the outflow 
has to adjust to the new situation. How fast this happens de- 
pends on the recombination time scale, which scales as ~ n"'. 
This method has been applied to a 280 ks XMM-Newton ob- 
servation of NGC 3783 with wildly different results. From the 
RGS data R >10 pc was inferred based on the lack of change 
in the deep O vii/O viii absorpti on edges and the Fe-M UTA ab- 
sorption complex dBehar et al.l 120031) . On the other hand, vari- 
ability in t he Fe xxv resonanc e Une at 6.7 keV detected in the 
EPIC data dReeves et al.l2004 impHes R <0.2 pc. From a 100 ks 
XMM-Newton observation of NGC 40 51, R ^ 0.001 pc wa s 
deduced for the highest ionisation gas (iKrongold et al.ll2007l) . 
However, the large amplitude variations (max/min =12) and 
the short time scales (down to 100 s) make these results rather 
model-dependent; an analysis of C handra LETGS data yiel ded 
distances in the range of 0.02-1 pc dSteenbr ugge et al.l l2009l) . 

Density-sensitive lines have given more robust results in the 
UV, and in at least one case (NGC 3783) they yield a dis- 
tance comparable to the one determined from densities based 
on recombin ation time scales. U sing metastable transitions in 
[Cm] ^1176, iGabel et al.1 d2005h find a distance of -25 pc for 
the absorbjing gas in one of the components in NGC 3783. 
They find a comparable distance based on the recombination 
time for Siiv. Again using the [Cm] metastable transitions, 
gas in one o f the components in NGC 4151 lies at a distance 
of <0.1 pc (Kraemer et al.' '2006"). Using metastable levels in 
[Sin and [Fen], MoeetaL (2009) find a distance for a low- 
ionisation outflow in the quasar SDSS J0838-I-2955 of ~ 3 kpc. 
Similarly, large distances for outflows have bee n found using 
similar density diagnostics in luminous quasars dHamann et al.l 



l200TllDunn et al.ll2010tlAoki et al .1120111 1 Arav et al.ll20'Tir 

In order to obtain a reliable density, we need both accurate 
ionic column densities, which yield the ionisation structure and 
total outflowing column density, and time variability on a suit- 
able time scale, from which the location of the outflow can be 
constrained through monitoring. With the current generation of 
X-ray telescopes, only a handful of AGN with outflows have 
high enough fluxes to yield suitable data for such an ambitious 
programme. Of these targets, Mrk 509 is the most promising due 
to the following three attributes. 1) It is among the brightest X- 
ray AGNs in the sky. 2) It has an excellent X-ray line str ucture 
dPounds et al.ll200lHSmith et al.ll2007HDetmers et aHomQ ). not 
as deep and blended as NGC 3783, but not as shallow and sparse 
as Mrk 279. Mrk 509 is a representative example of a moder- 
ate outflow. 3) Mrk 509 varies on a time scale of a few days 
making it ideal for a spectral timing campaign. Other targets 
like MCG -6-30-15, Mrk 766, or NGC 4051 with long XMM- 
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Newton exposure times vary too rapidly to do this: within the 
characteristic time scales of 100-1000 s, it is simply not possi- 
ble to get a high-quality grating spectrum. 

2.2. Kinematics of fhe outflow as traced in fhe UV 

Our campaign on Mrk 509 includes HST high-resolution ul- 
traviolet spectral observations that are simultaneously with 
Chandra's grating X-ray observations. The ultraviolet spectra 
provide a more detailed view of the kinematics of the outflow- 
ing gas. Prior high-resolution spectral observ ations with the Far 
Ultraviolet Spectroscopic Explorer (FUSE) dKriss et al.l l2000t) 
and the Space Telescope Imaging Spectrograph (STIS) on the 
Hubble Space Telescope (HST) (Kraemer et al. 2003) provide a 
good baseline for our observations, but these were both obtained 
nearly a decade before to our current campaign. The new UV 
spectra obtained with the Cosmic Origins Spectrograph (COS) 
on HST provide an updated view of the outflow, and enable us 
to examine long-term characteristics of the variability. 

Although the UV absorbers are not necessarily a per- 
fect tracer of the X-ray absorbing gas ("Kriss et al] 120001: 
iKraemer et al.ll2003l) . they provide insight into both lower ioni- 
sation absorbing components and portions of the high-ionisation 
outflow. Our COS observations yield accurate measurements of 
(or limits on) the column densities of C ii, C iv, N v, O i. Si ii - 
Si IV, S II, S III, Fe n, Fe iii, and H i Lyo-. Our limits or detections 
on the low-ionisation ions of Fe, Si, and S nicely complement 
the high-ionisation range covered by the RGS and the LETGS. 
Velocity-resolved measurements of the Li-like doublets of C iv 
and N v can be used to determine co lumn densities an d covering 
fractions of the outflowing gas (e.g.. lArav et al.ll2007l) . 

2.3. Abundances in AGN 

Studying the abundances, especially those of C, N, O, and Fe, 
of the gas in galaxies informs us about the ongoing enrichment 
processes through AGB star winds and supernova type la and 
II explosions. Therefore there has long been an interest in study- 
ing abundances with redshift to determine the star formation rate 
through history. Obvious sources for the study of the abundances 
are AGN, by taking advantage of their luminosity. Generally, the 
abundances in AGN have been determined from broad emission 
lines, however several assumptions make this method subject to 
systematic errors. Thanks to the improvement in atomic data, es- 
pecially in the dielectronic recombination rates for iron, the ex- 
cellent statistics in the combined RGS spectrum, and the simul- 
taneous optical and hard X-ray fluxes to constrain the spectral 
energy distribution (SED), we can determine reliable and accu- 
rate relative abundances from the narrow absorption lines of the 
outflow observed in Mrk 509. 

Since absolute abundances (i.e., metal abundances relative to 
hydrogen) cannot be measured directly from X-ray spectra, our 
UV observations of the Lyo' line provide a constraint on the total 
hydrogen column density. The first absolute abundance estimate 
for an outflow was derived by our group from the UV spectra 
of Mrk 279, where absorption lines of H, C, N, and O were 
carefully modelled (lArav et al.ll2007l) to find abundances rela- 
tive to solar for carbon (2.2+0.7), nitrogen, (3.5±1.1) and oxy- 
gen (1.6+0.8), which fully agree with the relative abundances 
of C, N, and O derived from the simultaneous X-ray spectra 
(ICostantini et al. 2007). The apparently enhanced N/H and N/O 
ratios in the Mrk 279 outflow may hint at strong contributions 
from stellar winds from massive stars and AGB stars. 



2.4. Broad emission iines 

Broad lines, especially at the energy of the Ovii triplet and 
the Oviii Lya line, have been recently detected in high- 
resolution X-ray spectra (Kaastra et al. 2002: Ogle et al. 20041: 
ISteenbrugge et al.ll2005t ICostantini et al.ll2007l: iLonginotti et all 
I20T0I). The width of the lines (about 1 A, 10000 kms"^) is com- 
parable to the width of the broad lines detected in the UV band 
and known to be produced in the so-called broad line region. 
In at least one case, the soft X-ray emission lines could be re- 
lated, v ia a physical model (t he locally optimall y emitting clouds 
mode l, [Baldwin et al.lll995h . to the UV lines (ICostantini et al.1 
l2007h . The contribution of the broad line region gas to the iron 
Kff line at 6.4 keV is not completely understood. For classical 
Seyfert 1 ga laxies this contribution seems modest (e.g. <20% 
in Mrk 279, ICostantini et al.l l2010l) . while it may account for 
most of the line emission in specific objects (e.g. NGC 7213, 
Bianchi et al. 2008). The quality of the data of Mrk 509 aflows us 
to significantly detect broad emission features around the O vii 
resonant line, the O viii, N vii Lyo- lines and the Ne ix triplet. The 
simultaneous observation of HST/COS and Chandra/LETGS al- 
low us to connect the broad lines detected in the UV with the 
soft energy lines and even the iron Ka line. 



2.5. Iron K complex 

Mrk 509 also shows a rich variety of emission and absorption 
components in the Fe-K band. In particular, EPIC data from pre- 
vious XMM-Newton observations shows evidence of a complex 
Fe-K emission line, with a narrow and neutral component pos- 
sibly produced far from the source, plus a broad, ionised, and 
variable componen t possibly originating in the accr etion disk 
(iPounds et albooUlPage et al.ll2003LlPonti et al.ll200 9). 

In addition, strong absorption features were found in the 
same data set at rest-frame energies 8-8.5 keV and 9.7 keV 
These were interpreted as being produced by H-like iron K and 
K-shell absorptions associated with an outflow with mildly rela- 
tivistic velocity of 0.14 - 0.2c. The lines were found to be vari- 
able in energy and marginal in intensity, implying that varia- 
tions in either the column density, geometry, and/or ionisation 
structure of the outflo w were maybe common in this source 
(Dadinaetal. 2005: Cappi et alJl2009l) . 

The above properties, combined with the source brightness 
(f2-i0kev ~ 2 - 5 X 10"'"* Wm-2) make Mrk 509 unique in 
attempting time-resolved spectral studies in the Fe-K band to 
follow the time evolution of both emission and absorption Fe-K 
features, and in disentangling the different physical components 
present in this source. 



2.6. Continuum emission of Seyfert galaxies 

The physical process at the origin of the X-ray emis- 
sion of Seyfert galaxies is generally believed to be thermal 
Comptonisation. In this process, the soft UV photons coming 
from the accretion disk are up-scattered into the X-ray domain 
by the hot thermal electrons of the corona. A thermal distribu- 
tion is generally preferred to a non-thermal one since the discov- 
ery, more than 15 years ago, of a high-en ergy cut-ofl" near 100 
keV in NGC 4151 b y OSSE and Sigma (iJourdainet al.l [19921: 
iMaisack et al.l [l993l) . Such a cut-off was not expected by the 
non-thermal models develope d at that time to e xplain the power- 
law-like X-ray spectrum (e.g. IZdziarski et al] 11990). The pres- 
ence of a high-energy cut-off has now been observed in a large 
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number of o bjects (iPerola et alJ l2002t iBeckmann et alJ l2005t 
lDadinall20Q8h . 

The cold (disk) and hot (corona) phases are expected to be 
radiatively Hnked: part of the cold emission, which gives birth to 
the UV bump, is produced by the reprocessing of part of the hot 
emission. Inversely the hot emission, at the origin of the broad- 
band X-ray continuum, is believed to be produced by Compton 
up-scattering of the soft photons, emitted by the cold phase, on 
the coronal energetic electrons. The system must then satisfy 
equilibrium energy balance equations, which depend on geom- 
etry an d on the ratio of direct he ating of the disk to that of the 
corona jHaardt & M araschi"1991). In the limiting case of a 'pas- 
sive' disk, the amplification of the Comptonisation process is 
only fixed by geometry. Therefore, if the corona is in energy 
balance, its temperature T^ and optical depth t must satisfy a re- 
lation that can be computed for different ge ometries of the disk 
-I- corona configuration (e.g.. lSvenssonl[l996l) . 

Realistic thermal Comptonisation spectra have been com- 
puted for more than two decades, and important effects like the 
anisotropy of the soft photon field have been precisely taken into 
account (Haardt 1993; Stern et al. 1995; Poutanen & Svensson 
Il996l) . These effects appear far from negligible. Noticeably, the 
spectral shape is significantly different for different disk-corona 
geometries but also for different viewing angles. Moreover, 
these effects underline the differences between realistic thermal 
Comptonisation spectra and the cut- off power-law appro xima- 
tion generally used to mimic them dPetrucci et alj|200lh . This 
may have important consequences in the study of superimposed 
spectral components like the soft X-ray excess and/or the the 
iron line and the reflection hump, which require a precise deter- 
mination of the underlying continuum. 

A well known characteristic of thermal Comptonisation 
spectra, however, is that they are strongly degenerate, i.e., signif- 
icantly different combinations of temperature and optical depth 
of the corona give the same power law slope in the 2-10 keV 
band. To break this degeneracy requires broad-band observa- 
tions from UV to hard X-rays, the X-ray/y-ray shape con- 
straining the high-energy cut-off, which is directly linked to 
the corona temperature. The ratio of the UV to soft X-ray flux 
constrains the optical depth better. Stronger and less ambigu- 
ous constraints on the nature of the coronal plasma can also 
be obtained from multiwavelength (from UV to hard X-rays/y- 
rays) variability studies, as they provide direct insight into the 
way the emitting particles are heated and cooled. Indeed, varia- 
tions in the X-ray spectral shape may be produced by intrinsic 
changes in the hot corona properties (e.g. changes of the heat- 
ing process efficiency) and/or by variations in the external en- 
vironment, such as changes in the soft photon flux (and con- 
sequently in the coronal cooling) produced by the cold phase 
(iMalzac & JourdainI I2OOOI: iPetrucci eTaD l2000l) . For example, 
thermal Comptonisation models, where hot and cold phases 
are in radiative equilibrium, predict that the X-ray spectrum 
of the sources should harde n when the energy of the high- 
energy cut-off increases (e.g. iHaardt et al.lil997.) . Indeed a cor- 
relati on between F and r, has been observed in di fferent ob- 
jects ( IPetrucci et al.ll2000l; IZdziarski & Grandill200ll) . thus giv- 
ing strong support to a thermal nature of the corona parti- 
cle distribution. The analysis of the one month, simultaneous 
lUE/RXTE monitoring campaign on NGC 7469, performed in 
1996 , is also in agreement with therma l Comptonisation emis- 
sion (Nandra et al 2000';'Petrucci et al."2004). 

The XMM-Newton/INTEGRAL monitoring of Mrk 509 pre- 
sented here provides an ideal data set to test Comptonisation 
models and to derive further constraints on the physical pa- 



rameters and geometry of the source, as well as on the precise 
shape of spectral components such as the soft-X-ray excess, the 
outflow, and the reflection component. It should also be noted 
that a good knowledge of the SED up to 100 keV is very im- 
portant input for the photo-ionisation modelling of the absorp- 
tion/emission features produced by these outflows. 

3. Mrk 509 

Mrk 509 was detected as a Seyfert 1 galaxy with a photographic 
magn itude mpg = 13 and a size of 10" ( Markarian & Lipove tskijI 
1 19731) . It has°a redshift of 0.034397 dHuchra et alj [1991 and 
broad (FWHM 4 500 kms ') and narrow optical emission lines 
(lMarkarianlll973h . Its optical luminosity is high for a Seyfert 1 
galaxy, and it puts Mrk 509 close to the limit between Seyfert 
galaxies and QSOs (Kopvlo vet al .1119741) . 

The host galaxy extends approximately east-west 



(iMagnitskaia & Saakian _ 1976h and has an axial ratio 



0.85 (iDahari & De RobertisI 119881; iKMiakos & Steineii 
Mrk 509 presently shows no nearby neighbours 



bla 

1990'). 

( Magnitskaia & SaakianI [l976h. the closes t spiral galaxy is at 
a pro jected distance of 0.3 Mpc dFuente s-Williams & Stocl3 
'1988|), but the asymmetric halo to the south (Mac KentViil99dh 
is direct evidence for strong interactions at some point in the 
recent past. 

In almost any wavelength band that opened its window, 
Mrk 509 has been one of the first AGN to be studied, thanks 
to its brightness, and several well-known astrophysicists today 
have investigated this object in earlier stages of their careers. 

In the radio b and, Mrk 509 was discovered at 3.9 cm in 1976 
dMingaliev et al.l [l978.) . Spatially resolved measurements at 6 
and 2 cm show an asyrametr i c structure of 1.4" x0.5" in PA 
124° dUlvestad & Wilsonlll984tlUnger et al.lll987h . correspond- 
ing to a linear size of about 1 kpc. This core is surrounded by 
more diffuse emission, visible at 20 cm, w ith a size comparable 
to the optical size (10 x 8", or 7 x 5 kpc ( Singh & WestergaardI 
Il992h . The radio emission seems to be preferentially aligned at 
PA -40° to -60°, coinciding with the preferred direction of the 
optica l polarisation ( Singh & Westergaard 1992; Martin et aT| 
119831) . 

In the infrared the first observati ons date back to 1975 
dStein & WeedmanI 119761; lAUenl Il976l) . The spectrum in the 
l-lO/zm band is approximately a pow er law ( Riekd 119781 ; 
iGlass et alJ Il982t iMcAlarv et all Il983h . Contrary to several 
other Seyferts, Mrk 509 shows no evidence of strong dust fea - 
tures in infrared spectra dRoche et al.lll984t lMoorwoodlll986l) . 
The variability in the mid-IR is much weaker than at shorter 
wavelengths, with a possible delay of two to three months (Glass] 
2004). This suggests that the bulk of the IR emission originates 
far from the nucleus. 

Quantitative measurements of the optical broad and narrow 
emission lines were first p resented by Osterbrock (1977) and 
Ide Bruvn & Sargena dl978h . Since then, dozens of papers have 
appeared with refined measurements and models for the line pro- 
files and intensity ratios. The size of the broad emission line re- 
gion in the Wfi line has been estimated as 80 light days (0.07 pc) 
based on reverberation studies (Peterson et al. 1998). 

Perhaps the most detailed study of interest for our project 
has be en the mapping of the [Oiii] /15007 line bv .Phillips et al.l 
dl983h . They show that the line has two components. There 
is a rotating, low-ionisation gas disk, coinciding with hot stel- 
lar components, which is photo-ionised by the UV radiation of 
young hot stars. The rotation axis has approximately the same di- 
rection as the radio emission and the preferred polarisation angle 
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of the optical emission mentioned before. The second compo- 
nent is expanding, high-ionisation gas, which is photo-ionised by 
the nucleus. The line profile of that component is consistent with 
outflowing gas up to velocities of -800 km s"' . It extends out to 
about 5 kpc from the nucleus. Outflowing gas is seen at opposite 
sides of the nucleus, indicating that the outflow is approximately 
face-on. A similar spatial extent is seen in a few other narrow 
emission lines (Winge et al. 2000). Approximately the same re- 
gion also shows broad Hfi line emission, most likely caused by 
scattered nuclear light in the narrow line region (iMediavilla et al.l 
Il998l) . X-rays from the nucleus travel through these media, and 
the lowewst ionised X-ray absorbers may show the imprint of 
this material on the spectrum. 

The first UV spectra were presented in 1980, with measure- 
ments of the Lya Une flux using the lUE satelhte dWu et al.l 
Il980l) . Sig nificantly better spect ra we re subsequently obta ined 
with HST ("Crenshaw e t all 19951 [19991 iKraemer et al.l2003l) and 
FUSE (Kriss et al...2000l) ! 

In X-rays, the first detection was made with A riel V in 
1974-1976 in the 2-10 keV band ("Cooke et al."1978'). Soon the 
first (power-law) spectra were fitted (Mushotzkv et al. 1980), 
and the source w as detected at higher energies, beyond 40 keV 
(lDiletal.1119811) and u p to about 200 keV (.Pietschet al...l981; 
iRothschild et al.|[l983h . It was one of the first sour ces where a 
soft X -ray excess was discovered with HEA01-A2 ('Singh et al.' 
Il985h. the first iron line detection dated to 1987 with EXOSAT 
dMorini e t al. 198 71), and the r eflection component was revealed 
by Ginga (Pounds etal.ll 1994 . 

Variability was discover ed first in the optical band 
(^Magnits kaia & Saakianlll976l) . They found variations of about 
a magnitude on a time scale of several months, but a compar- 
ison with observations 21 years before showed a similar flux 
level. This is one of the desired properties for the purpose of 
our campaign: significant variability on suitable time scales, 
with a limited range of flux levels. Since then, variability has 
been found in all energy b ands from r adio (Mingaliev et al. 
1 19781). infrared co ntinuum (iGlassI l2004l) , optical polarisation 
("Martin et al."1983), broad emissi on Unes dPeterson et alj l 19821 
[1984; Rosenblatt et al. 1992), UV (IChapma n et al. 1985), and X- 
rays (iDoweret al.lll980l) . 

There have been many papers that derive the mass of the cen- 
tral SMBH using a variety of methods. Derived numbers range 
from 1.4 X 10^ M© in one of the oldest determinations, based 
on the co rrelation of the H/S intensity with the full-width at zero 
intensity (lLiulll983h to 1.43 + 0.12 x 10*^ Mq for a modern es- 
timate based on reve rberation mapping of optical broad lines 
dPeterson et al.ll2004 . 

4. Observations 

We have obtained data from seven different observatories during 
our campaign. At the core of our programme are ten observa- 
tions of approximately 60 ks each, with XMM-Newton, spaced 
by four days. We used the data of all instruments on XMM- 
Newton: Reflection Grating Spectrometer (RGS), EPIC (pn and 
MOS), and Optical Monitor (OM). For this last instrument, all 
filters except for the white light filter were used, and we also ob- 
tained a spectrum with the optical grism for each observation. 
Observations with the UV grism were not allowed for opera- 
tional reasons. 

Simultaneously with our ten XMM-Newton observations, 
we obtained data with INTEGRAL to observe the hard 
X-rays. Three weeks after the end of the XMM-Newton 
monitoring, we did simultaneous observations with Chandra 
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Fig. 1. Timeline of our monitoring campaign of Mrk 509. The 
first observation, with Swift, started on September 4, 2009, and 
the last observation, with Chandra, ended on December 13, 
2009. 

(Low Energy Transmission Grating (LETGS) with the HRC- 
S camera), together with the Hubble Space Telescope (HST) 
Cosmic Origins Spectrograph (COS). Due to operational con- 
straints, these observations could not coincide with the XMM- 
Newton/INTEGRAL monitoring. 

Before our XMM-Newton/INTEGRAL observations, and in 
between the XMM-Newton and Chandra/HST observations, we 
monitored Mrk 509 with Swift, using both its X-ray telescope 
(XRT) and the Ultra Violet and Optical Telescope (UVOT). This 
was in order to obtain a continuous monitoring and to allow us 
to study the flux before the start of our campaign. This is impor- 
tant because absorption components may respond with a delay 
to continuum variations. 

Finally, throughout our high-energy campaign, we obtained 
optical spectroscopy and photometry. We have five observations 
with the 4.2 m William Herschel Telescope (WHT) at La Palma, 
photometry with the ACAM camera using the Sloan g, r, i and Z 
filters, and low-resolution spectroscopy with the VPH disperser 
Furthermore, we obtained nine observations with the 1.3 m 
PAIRITEL (Peters Automated IR Imaging Telescope), with pho- 
tometry in the J, H, and K bands. Table[T]gives some details on 
our observations, and in Fig. [T] we show a graphical overview of 
the timeline of our campaign. 

5. Light curve 

The light curves in the X-ray band (0.3-1.0 keV, 2-10 keV and 
20-60 keV) and in UV (2310 A) are shown in Fig. |2] More 
details and the interpretation of these l ight curves are given 
in subsequent pape rs of our series (e.g. iMehdipour et alJfeOl lb 
iPetrucci et al.ll201 II) . The flux was at a typical flux level, com- 
pared with archival observations. Right in the middle of our ob- 
serving campaign, the source showed a significant brightening 
by ~ 50%, in particular in the soft X-ray band, and with a smaller 
amplitude in the UV and hard X-ray bands. This allowed us to 
study the response of the absorption components to this bright- 
ening, and thanks to our broad-band measurements we can con- 
strain the emission mechanisms for the continuum components. 



6. Spectral energy distribution 

The broad-band spectral energy distribution (SED) is essen- 
tial for obtaining the ionisation balance needed for the photo- 
ionisation modelling of the outflow. The multiwavelength nature 
of our campaign gave us almost simultaneous coverage of the 
total spectrum, which allowed us to constrain the exact shape of 
the SED to high accuracy. 

We need two different SEDs for our campaign. The first de- 
scribes the average SED during the XMM-Newton part of the 
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Table 1. The observation log of our Mrk 509 campaign. 









Start time (UTC) 


Span 








Start time (UTC) 


Span 


Telescope 


Obs. 


ID 


yyyy-mm-dd@hh:mm 


(ks) 


Telescope 


Obs. 


ID 


yyyy-mm-dd@hh:mm 


(ks) 


Chandra 


1 


11387 


2009-12-10@04:54 


131.4 


HST/COS 


1 


IbdhOlOlO 


2009-12-10®02:48 


2.0 


Chandra 


2 


11388 


2009-12-12@22:33 


48.6 


HST/COS 
HST/COS 


2 
3 


lbdh01020 
lbdh01030 


2009-12-10®04:07 
2009-12-10®05:42 


2.7 


XMM-Newton 


1 


0601390201 


2009-10-15@06:19 


57.1 


5.5 


XMM-Newton 


2 


0601390301 


2009-10-19® 15:20 


59.9 


HST/COS 


4 


lbdh01040 


2009-12-10®08:54 


2.7 


XMM-Newton 


3 


0601390401 


2009-10-23@05:41 


60.5 


HST/COS 


5 


lbdh02010 


2009-12-1 1®02:46 


2.0 


XMM-Newton 


4 


0601390501 


2009-10-29@06:55 


60.5 


HST/COS 


6 


lbdh02020 


2009-12-1 1®04:05 


2.7 


XMM-Newton 


5 


0601390601 


2009-1 1-02@02:46 


62.4 


HST/COS 


7 


lbdh02030 


2009-12-1 1®05:41 


2.7 


XMM-Newton 


6 

7 


0601390701 
0601390801 


2009-1 1-06@07:00 
2009-ll-10@08:42 


62.7 
60.5 


HST/COS 


8 


lbdh02040 


2009-12-ll®07:16 


5.5 


XMM-Newton 


WHT 


1 


1372032-154 


2009-10-04®22:50 


1.6 


XMM-Newton 


8 


0601390901 


2009-1 1-14@08:27 


60.5 


WHT 


2 


1372732-166 


2009-10-07®20:10 


3.6 


XMM-Newton 


9 


0601391001 


2009-1 1-18@02:08 


65.1 


WHT 


3 


1378102-127 


2009-1 1-01 ®20:55 


2.4 


XMM-Newton 


10 


0601391101 


2009-1 1-20@07:40 


62.4 


WHT 
WHT 


4 

5 


1378849-900 
1382682-702 


2009-1 1-06® 19:30 
2009-1 1-27® 19:42 


5.4 


INTEGRAL 


1 
2 


07200160001 
07200160002 


2009-10-14@21:55 
2009-10-19® 19:38 


127.7 
120.1 


1.5 


INTEGRAL 


Swift 


1 


00035469005 


2009-09-04® 13: 18 


0.9 


INTEGRAL 


3 


07200160003 


2009-10-22@ 19:29 


124.3 


Swift 


2 


00035469006 


2009-09-08®02:30 


1.5 


INTEGRAL 


4 


07200160004 


2009-10-28@ 19:08 


110.7 


Swift 


3 


00035469007 


2009-09-12® 18:52 


0.5 


INTEGRAL 


5 


07200160005 


2009-ll-01@01:07 


122.3 


Swift 


4 


00035469008 


2009-09-16®20:36 


1.1 


INTEGRAL 


6 


07200160006 


2009-11-05® 15:00 


63.0 


Swift 


5 


00035469009 


2009-09-20® 07:51 


1.1 


INTEGRAL 


7 


07200160011 


2009-11-06® 18:32 


59.0 


Swift 


6 


00035469010 


2009-09-24® 09:51 


1.3 


INTEGRAL 


8 


07200160007 


2009-ll-09@18:18 


120.7 


Swift 


7 


00035469011 


2009-10-02® 18:45 


1.0 


INTEGRAL 


9 


07200160008 


2009-ll-13@20:13 


122.3 


Swift 


8 


00035469012 


2009-10-05®23:59 


1.0 


INTEGRAL 


10 


07200160009 


2009-11-17® 14:48 


61.2 


Swift 


9 


00035469013 


2009-10-10®21:37 


1.1 


INTEGRAL 


11 


07200160010 


2009-11-19® 11:23 


110.7 


Swift 
Swift 


10 
11 


00035469014 
00035469015 


2009-10-14® 10:09 
2009-10-18®05:56 


1.4 


PAIRITEL 


1 


166.2 


2009-09- 15® 04:45 


0.8 


1.0 


PAIRITEL 


2 


166.2 


2009-09-16®04:35 


0.8 


Swift 


12 


00035469016 


2009-ll-20®07:18 


1.0 


PAIRITEL 


3 


166.3 


2009-09-20® 04:51 


0.8 


Swift 


13 


00035469017 


2009-1 1-24® 23:39 


1.0 


PAIRITEL 


4 


166.4 


2009-09-24®03:55 


0.8 


Swift 


14 


00035469018 


2009-ll-28®19:13 


1.2 


PAIRITEL 


5 


166.5 


2009-09-28®03:24 


0.8 


Swift 


15 


00035469019 


2009- 12-02® 00: 18 


1.3 


PAIRITEL 


6 


166.7 


2009-10-22®02:55 


0.8 


Swift 


16 


00035469020 


2009-12-06®07:07 


1.2 


PAIRITEL 


7 


166.8 


2009-10-27®04:06 


0.8 


Swift 


17 


00035469021 


2009-12-08®02:30 


1.1 


PAIRITEL 


8 


166.9 


2009-10-31 ®03:21 


0.8 


Swift 


18 


00035469022 


2009-12-10® 12:36 


1.2 


PAIRITEL 


9 


166.10 


2009-ll-17®03:19 


0.8 


Swift 


19 


00035469023 


2009-12-12®20:30 


1.0 



campaign, which is needed for the analysis of the stacked RGS 
spectrum. The second one is needed for the Chandra LETGS and 
HST COS part of the campaign, when the source had a lower flux 
level. Table |2] shows the flux points of every instrument used to 
create the SED shown in Fig. [3] The statistical uncertainties on 
the flux are much lower than the reported values (except in the 
case of INTEGRAL), therefore we list the systematic uncertain- 
ties due to calibration uncertainties. 



6.1. X-ray SED 



We u sed the best fit RGS model (model 2 of iDetmers et al.l 
12011 ) for the spectrum between 7-38 A. We correct for both 
the intrinsic absorption (the ionised outflow) and the Galactic 
ISM absorption. For the flux between 1.2 A (10 keV) and 7 A 
(1.77 keV), we used the EPIC-pn data. Above 10 keV we used 
the INTEGRAL data, which were obtained simultaneous with 
the XMM-Newton observations. This gives us the continuum 
flux up to ~ 200 ke V. B eyond that, the SED is extrapolated using 
the model described in lPetrucci et al.l (1201 ll) . 

We compared the Chandra LETGS flux in two different en- 
ergy bands with the non-contemporaneous RGS flux to obtain 



the flux variations in the soft (20 - 35 A) and hard (7 - 10 A) X- 
ray bands. The RGS flux in these bands is 66% and 30% higher 
than the LETGS flux, respectively. The soft X-ray flux for the 
RGS observation for A > 31 A is obtained by scaling the LETGS 
continuum by a factor of 1 .66. As we have no information on the 
X-ray flux above 10 keV during the LETGS observation, and the 
uncertainty on the LETGS spectrum for A < 1 A increases, we 
estimate the flux for the Chandra observation for /t < 7 A by 
dividing our model for XMM-Newton/INTEGRAL by a factor 
of 1.30. 



6.2. EUV extrapolation 

The EUV spectrum produces most of the ionising flux; however, 
it is also the most uncertain part of the SED, as there are no data 
between 50 and 912 A. LETGS formally measures up to 175 A, 
but due to the Galactic absorption, the flux is low and also the 
modelling of the higher spectral orders becomes more uncertain 
at longer wavelengths. We therefore have to interpolate our SED 
between the soft X-rays and the UV. There are several options 
for doing this. Our baseline model is a power-law interpolation 
between 50-943 A (the last FUSE data point). For the XMM- 
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Table 2. Continuum fluxes corrected for Galactic and intrinsic 
absorption during the XMM-Newton and Chandra parts of the 
campaign. 
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MJD - 55075 (days) 

Fig. 2. Top panel: UV light curve of Mrk 509 as obtained with 
Swift/UVOT, XMM-Newton/OM, and HST/COS. The data have 
been corrected for extinction. Bottom panel: X-ray light curve 
as obtained by Swift/XRT, XMM-Newton/pn, Chandra/LETGS, 
and INTEGRAL/ISGRI. These data have also been corrected for 
Galactic absorption. For the pn data, the squares correspond to 
the 0.3-1.0 keV band, the stars to the 2-10 keV band. Note that 
10'2 JyHz equals lO"'"* Wm^^ (IQ-" in c.g.s. units). 
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Fig. 3. The SED used of Mrk 509 for the time-averaged XMM- 
Newton and Chandra observations. Data points from Table|2]are 
indicated with triangles and squares. 



Newton epoch this power law has a photon index of 2.39, for the 
Chandra epoch 2.53. Alternatively, we consider some kind of 
big blue bump, by extrapolating the RGS spectrum with a power 
law, using its slope near 30 A (photon index 2.58) up to the point 
(193 A) where it matches the extrapolation from the UV with a 
photon index of -2. The different extrapolations between the X- 
ray and UV data are shown in Fig. [3] 

6.3. UVSED 

The UV part of the SED was derived from archival FUSE data 
and from the HST COS data of our campaign. As FUSE and 
COS share a common wavelength band, the FUSE archival 



Frequency 


E/A 


vF„ " 


vF." 


Unc.'- 


Instr.'' 


(Hz) 








(%) 




3.91x10''' 


162 keV 


2.67 


2.05 


50 


INTEGRAL 


2.79x10''' 


115 keV 


3.08 


2.37 


30 


INTEGRAL 


2.20x10''' 


91keV 


3.39 


2.60 


10 


INTEGRAL 


1.57x10''' 


65keV 


3.72 


2.86 


10 


INTEGRAL 


1.05x10''' 


43keV 


4.16 


3.20 


10 


INTEGRAL 


7.30x10"* 


30keV 


4.37 


3.36 


10 


INTEGRAL 


4.84x10"* 


20keV 


4.50 


3.46 


10 


INTEGRAL 


2.00x10"* 


1.5 A 


3.35 


2.58 


5 


pn 


1.20x10"* 


2.5 A 


2.95 


2.27 


5 


pn 


6.00x10'^ 


5.0 A 


2.58 


1.98 


5 


pn/RGS 


3.00x10'^ 


10 A 


3.00 


2.21 


5 


RGS 


2.00x10'^ 


15 A 


3.61 


2.51 


5 


RGS 


1.50x10'^ 


20 A 


4.33 


2.78 


5 


RGS 


1.20x10'^ 


25 A 


5.25 


3.24 


5 


RGS 


1.00x10'^ 


30 A 


5.72 


3.46 


5 


RGS 


8.57x10'^ 


35 A 


6.46 


3.89 


5 


RGS 


7.50x10'* 


40 A 


6.83 


4.12 


5 


LETGS 


6.67x10'* 


45 A 


7.08 


4.27 


5 


LETGS 


6.00x10'* 


50 A 


7.21 


4.35 


5 


LETGS 


3.18x10'^ 


943 A 


22.6 


20.8 


10 


FUSE 


3.13x10'^ 


960 A 


20.6 


19.0 


10 


FUSE 


3.05x10'^ 


984 A 


20.6 


18.8 


10 


FUSE 


3.02x10'^ 


993 A 


19.0 


17.5 


10 


FUSE 


2.55x10'^ 


1175 A 


22.1 


20.3 


5 


FUSE 


2.12x10'^ 


1415 A 


23.3 


21.5 


5 


COS 


1.71x10'^ 


1750 A 


20.8 


19.1 


5 


COS 


1.42x10'^ 


2120 A 


18.2 


17.5 


3 


OM 


1.30x10'^ 


2310 A 


17.6 


17.1 


3 


OM 


1.03x10'^ 


2910 A 


13.9 


13.5 


3 


OM 


8.72x10'" 


3440 A 


12.4 


12.1 


3 


OM 


6.67x10'" 


4500 A 


9.24 


8.98 


3 


OM 


5.52x10'" 


5430 A 


7.65 


7.46 


4 


OM 


2.44x10'" 


1.23 yum 


3.17 


10 


V 


1.81x10'" 


1.66 yum 


4.01 


10 


W 
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5.30 


10 
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4.56 


15 
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15 
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For the time-averaged XMM-Newton spectrum, in units of 

10-'" Wm-2 (or 5.034 photons m-^ s-' A-') 

For the time-averaged Chandra spectrum, in units of 10'" Wm-^ 

(or 5.034 photons m-- s-' A-') 

Estimated systematic and statistical uncertainty 

Instrument used to derive the flux 

Values obtained from the literature for the J, H, K, IRAS and radio 

bands; see text for references. 



fluxes were scaled to the COS flux level, assuming that the spec- 
tral shape remained the same. 

The full description of the COS data reduction can be found 
in lKrissetaH (l201lh . In short, the COS data were flat-field cor- 
rected and additional wavelength calibration performed. Time- 
dependent sensitivity corrections were applied to the COS data, 
resulting in an absolute flux accuracy of 5%. As the COS data 
were taken during the LETGS observations, we needed to adjust 
them to the higher flux levels during the XMM-Newton observa- 
tions. The OM data at somewhat longer wavelengths (see below) 



8 



J.S. Kaastra et al.: Multiwavelength campaign on Mrk 509 



show on average 6% higher flux during flie XMM-Newton ob- 
servations compared to the LETGS observations. Using the near- 
simultaneous FUS E/HST/optical spectrum of Mrk 509 from 
IShang et aP (l2005i) normaUsed to the OM flux, we find that in- 
creasing the COS flux by 10% is appropriate if assuming a con- 
stant shape for the UV spectrum. 

6.4. Optical SED 

The optical part of the SED was obtained from the OM data. The 
data were corrected for interstellar absorpt ion an d de-reddened 
using the reddening curve of Cardelli et al.| (Il989l) . including the 
near-UV update by lO'DonnelL iil994,) . The total colour excess 
E{B - V) is 0.057 mag, and Ry = A{V)/EiB - V) was fixed at 
3.1. The total neutral hydrogen column densit y in the d i rectio n 
of Mrk 509 is 4.44 x 10^^ m'^, as given by Mur phy et al](ll995h 
The host gala xy correction is ba sed on the results of Ben tz et aP 
(l2009i) and iKinnevetalJ (Il996l) . For Mrk 509 the host-galaxy 
contribution at 2231 A is 7 x 10"^° Wm"^ A ', which is neg- 
ligible compared to the AGN flux of 8.45 x 10"^^ Wm"^ A"'. 
At longer wavelengths the host galaxy contribution increases, 
to 3.2 X 10"'** Wm-^A-' at 5500 A, where the AGN flux is 
1.7 X 10 '^ Wm"^ A '. Fo r a more thorough desc ription of the 
optical data reduction, see Mehdipour et al.l (1201 ih . The optical 
fluxes (from Swift UVOT) at the time of die LETGS observa- 
tions have to be adjusted for the lower flux level of the source. As 
XMM-Newton observation 2 has a similar flux to that of Swift 
UVOT observation 18, which is close in time to the LETGS ob- 
servations, we assume that the optical flux during the LETGS 
observations is the same as that of XMM-Newton OM observa- 
tion 2, which is 6% smaller than the average optical flux during 
the XMM observations. 



6.5. Infrared and radio SED 

At IR wavelengths shortward of 1 //m, the spectrum shows a 
small upturn due to emis sion from the torus (see, e.g., the spec- 
tra by iLandt et al.ll201l1) . The infrared fluxes that we use here 
are based on photometry with host-galaxy subtraction in the J, 
H, and K bands. W e use the average from the o bservations of 
iDanese et aTl (Il992h and lKotilainenetal.1 ( |1992|) . which are in 
good agreement with ea ch other For 100 //m, we use the IRAS 
flux (lMoshu-etal.lll990h . Between the K-band and the 100 yum 
flux points, the SED shows an almost consta nt value in term s 
of vFy. In the 6-35 yum band Spitzer data exist (W u et al.ll2009h . 
with ~ 30% higher flux compared to our adopted SED. However, 
these data, taken with a relatively large aperture of 10", are not 
corrected for the stellar contribution. As our results are not very 
sensitive to the details of the IR spectrum, we can safely use our 
simplified SED in the mid-lR band as shown in Fig. [3] 

Beyond 100//m, the flux drops rapidly; following usual prac- 
tice in this region, we extrapolate the flux from 100 fim to lower 
frequencies with an energy index of -2.5 down to 1 cm. This 
gives a flux of onl y about a factor of tw o above the observed 
value at 14.9 GHz (iBarvainis et al.|[l996h . Finally, we also con- 
sider an SED without IR emission, because it is unclear how 
much of the IR emission is seen by the outflow (see next sec- 
tion). 



6.6. Effect of the SED on the ionisation balance 

Because a different SED can have a strong eff'ect on the ionisa- 
tion balance of the outflow, we investigated the effects of differ- 
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Fig. 4. The four different cooling curves for the different as- 
sumed SEDs. 



ent SEDs, focusing on the uncertainties in the EUV and infrared 
parts of the spectrum. In the EUV band we have no measure- 
ments, so we have to rely on interpolation. Depending on the 
location of the outflow relative to the dusty torus expected to sur- 
round the AGN, it may receive either all or only a small fraction 
of the IR flux emitted by the torus. To investigate the effects of 
different assumptions about this, we compare four SEDs (Fig.O: 

1. Our standard SED (Table|2] column 3 and Fig.|2] solid line), 
which uses a simple power-law interpolation between the 
soft X-ray and UV bands; 

2. A SED with a stronger EUV flux, using the broken power- 
law approximation with a break at 193 A mentioned in 
Sect. 16.21 to mimic a stronger EUV flux (EUV bump in 
Fig.0; 

3. Same as model 1, but with the infrared flux essentially set to 
zero, for an absorber that does not receive emission from the 
torus (No IR in Fig.[3]i; 

4. the SED for the Chandra observation (Table |2] column 4), to 
see the effects of time variability. 

The ionisat ion balance c a lculat ions were performed us- 
ing Cloud>13 (iFerland et al.l Il998h version C08.00, with 
Lodders et al. (l2009t) abundances. The results are shown in 
Fig. |4] As can be seen, the differences between the four cool- 
ing curves are small. The case without infrared flux (model 3) 
starts deviating from the default case above 10^ K, because of 
enhanced cooling by inverse Compton scattering of the infrared 
photons. The case with a stronger EUV flux (model 2) has a sim- 
ilar shape to model 1, but has on average log H higher by 0.08 
due to the enhanced ionising EUV flux. This results in a small 
shift towards the right of the figure. Finally, the lower flux during 
the Chandra observations leads to a very similar cooling curve, 
because the peak of the UV spectrum is very similar 

We are therefore confident that our derived SED is an accu- 
rate description of the true source continuum of Mrk 509 and 
that the photo-ionisa tion results obtained from the analysis of 
iDetmers et al.l (12011) are not seriously affected by the uncertain- 
ties in the assumed SED. 

As an additional test we checked how metallicity influences 
the cooling curve. We did this because we do not know a priori 
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Fig. 5. The two different cooling curves for solar (solid line) and 
3 times solar metallicity (dashed line). 



the metallicity of the outflow, which can have a significant effbct 
on th e shape of the cooling curve (see e.g. IChakravortv et al.l 
[2009'). We adopted 3 times t he sola r abundances, cf. typical 
abundances found in Mrk 279 jArav e t al. 2007). Figure|5]shows 
the results compared to the solar metallicity run. There are slight 
differences in the shape of the two cooling curves, with the 
higher metallicity curve having bigger unstable branches (where 
the slope is negative) due to enhanced line cooling. However, 
the resulting photo-ionisation model will give no significant dif- 
ferences in derived parameters, except for the hydrogen column 
densities, which are a factor of three lower owing to the three 
times higher metal abundances. This again shows that the ioni- 
sation structure of the outflow in Mrk 509 is very stable to the un- 
certainties that are present in either the assumed SED or metal- 
licity. 

7. Summary 

In this paper we have introduced our multiwavelength study 
of the Seyfert 1 galaxy Mrk 509, which used XMM-Newton, 
INTEGRAL, Chandra, HST, Swift, WHT and PAIRITEL. Our 
observations spanned 100 days of monitoring from September 
2009 to December, covering band-passes from 2 //m to 200 keV. 
The core of our programme consisted of ten simultaneous ob- 
servations with XMM-Newton and INTEGRAL, followed by a 
long Chandra LETGS spectrum obtained simultaneously with an 
HST/COS far-UV spectrum. The high-resolution spectroscopy 
combined with the time variability monitoring enabled us to dis- 
entangle the different absorption and emission components in 
Mrk 509. 

Using our data we produced a continuous light curve in the 
X-ray and UV bands and a comprehensive SED. Thanks to our 
broad-band coverage, we deduced that uncertainties on the SED 
have little effect on the photo-ionisation equilibrium that applies 
to our subsequent models of the ionised outflow. Our light curve 
shows a 60% flux increase in the soft X-ray band correlated with 
an enhancement of the UV flux. 

A series of subsequent papers will elaborate on these re- 
sults and others from our campaign. The stacked XMM-Newton 
RGS spectrum is presented in Kaastra et al. ( 2011, paper II). 
The time-averaged ionisation and velocity struct ure of the out- 
flow deduced from this spectrum is presented by iDetmers et al.l 



(1201 ll paper III), the broad-ba nd continuum varia bility by 
'Mehdipour et all (1201 ll paper IV). 'Ebrero et al.' ("2011' paper V) 
describes the Chandra LETGS data and Kriss et al. (2011, pa- 
per VI) presents the analysis of the HST/COS data. More papers 
are in preparation. 
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